No-passing Rule in the Ground State Evolution of the Random-Field Ising Model 
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We exactly prove the no-passing rule in the ground state evolution of the random-field Ising model 
(RFIM) with monotonically varying external field. In particular, we show that the application of 
the no-passing rule can speed up the calculation of the zero-temperature equilibrium M{H) curve 
dramatically. 
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I. INTRODUCTION 

The no-passing rule was first introduced by Middleton 
in the study of sliding charge-density waves (CDW's) Q. 
The CDW problem belongs to the more general class of 
driven elastic manifolds in random media. If one de- 
fines a simple one-dimensional order parameter within 
the model, then a natural partial ordering of the con- 
figurations can be defined. In the simple CDW model 
considered by Middleton, the CDW configuration {(pi{t)} 
describes the CDW distortions at N lattice sites indexed 
by i, with ipi{t) real phase variables and t the time. The 
equation of motion for an overdamped CDW is governed 
by the Langevin dynamics 



= AV. - Va^r) + fit) 



(1) 



Here, the A'^ipi term represents the simple elastic in- 
teractions. Vi'{(pi) is the pinning force at site i due 
to the 2tt periodic pinning potential V{(pi). And f{t) 
stands for the external driving force. Then one can de- 
fine the natural partial ordering of two configurations: 
= {iff} > = if > foj. g^^i^ 

site i of the system. The no-passing rule states that 
given a driving force / if initially C^(0) > C^(0), then 
C^(i) > C^it) for all t > 0, i.e. the "greater" (C°) is 
never passed by the "lesser" (C^). As stressed by Mid- 
dleton, this rule relies crucially on the elastic potential 
being convex. In other words, the elastic potential tends 
to decrease the separation of nearest-neighbor (p's. More 
recently, Krauth et al found a similar no-passing rule 
in the study of driven elastic strings in disordered me- 
dia [1, Q. Obviously, this is the same general problem. 
Again, the rule is crucially dependent on the fact that 
the elastic potential is convex. 

The no-passing rule can be used to prove many use- 
ful properties, such as the asymptotic uniqueness of the 
sliding state for CDW's ^ and the intriguing memory 
effects Q- In fact, just after its introduction by Middle- 
ton, the no-passing rule was used in the non-equilibrium 
zero-temperature random-field Ising model (RFIM) by 
Sethna et al to prove the return point memory effect [J| . 
The RFIM is obtained by adding a random field hi at 
each site of the Ising model 



H - ^ J SiSj - (g -I- hi) Si 



(2) 



<ij> 



The distribution of hi values is usually taken to be Gaus- 
sian, with standard deviation R and mean 0. J is the 
nearest-neighbor ferromagnetic coupling strength and H 
is the uniform external field. In this case, the natural 
partial ordering of two configurations can be defined sim- 
ilarly as in the CDW case. A difference is that Si = ±1 
while ipi is real. The no-passing rule states: Let a sys- 
tem C^(t) be evolved under the fields H'^{t) and sim- 
ilarly C^{t) evolved under H^{t). Suppose the fields 
H'^{t) > H^{t) and the initial configurations satisfy 
C°(0) > C^(0), then C^(t) > C^[t) at all times later 
t > 0, i.e. the partial ordering will be preserved by the 
dynamics. With a local metastable single-spin-fiip dy- 
namics, i.e. a spin flips when its effective local field 



hf = J^s,+h, + H 



(3) 



changes sign, the proof of the no-passing rule is straight- 
forward ;4]. Even with a two-spin-flip dynamics, it has 
been shown by Vives et al that the no-passing rule is 
still true at zero-temperature Q . Note that for the mag- 
netization process, the no-passing rule is equivalent to 
the fact that the flipped spins can never flip back as H is 
swept monotonically. Again, this rule is not uncondition- 
ally true. It relies crucially on the nearest-neighbor inter- 
action being ferromagnetic (J > 0). Just like the convex 
elastic potential, the ferromagnetic interaction also tends 
to decrease the separation of nearest-neighbor degrees of 
freedom, i.e. it tends to align the spins. 

Recently, in the study of the equilibrium zero- 
temperature RFIM, Vives et al conjectured that when 
the external field H is swept from — oo to oo, flipped spins 
cannot flip back in the equilibrium M{H) curve @. In 
other words, the no-passing rule is valid even for the zero- 
temperature equilibrium dynamics, i.e. the evolution of 
the ground state (GS). Vives et al further conjectured 
that this rule can be used to speed up the calculation 
of the equilibrium M[H) curve since flipped spins at a 
lower field can be removed from the GS calculation for 
all higher fields. Unfortunately, this simple but powerful 
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rule has not been proven so far for the equihbrium RFIM. 
This is the main motivation of our work. 

This paper is organized as follows. In Sec. II, we give 
a short introduction to the calculation of the equilibrium 
M{H) curve of the zero-temperature RFIM. In Sec. Ill, 
we work out some basic steps for the proof of the equi- 
librium no-passing rule. In Sec. IV, we present the proof. 
In Sec.V, we show the direct application of this rule to 
the calculation of the equilibrium M{H) curve. Finally, 
in Sec. VI we discuss its validity in other systems. 



II. THE EQUILIBRIUM M{H) CURVE 

To calculate the equilibrium M (H) curve of the zero- 
temperature RFIM, we first need to calculate the exact 
GS in the RFIM at an arbitrary applied external field 
H. This is the basic step of calculating the equilibrium 
M{H) curve, i.e. the GS evolution for varying H . For- 
tunately, there is a well-known mapping of the RFIM 
GS problem to a min-cut/max-flow problem in combi- 
natorial optimization. The mapping and the so-called 
push-relabel algorithm for the min-cut/max-flow prob- 
lem has been well described in the literatures 0, [1] ■ For 
RFIM, the run time of the push-relabel algorithm scales 
as 0{N^/^) with N the system size 

The equilibrium M{H) curve can be simulated with 
the method reported in Ref. [111. It is essentially based 
on the fact that the GS energy E{Si, H) is convex up in 
H, which allows for estimates of the fields H where the 
magnetization jumps (called "avalanches" occur). This 
algorithm finds steps by narrowing down ranges where 
the magnetization jumps with an efficient linear interpo- 
lation scheme. An illustration of the algorithm is shown 
in Fig. [TJ An example of the calculated equilibrium 
M{H) curve is shown in Fig El The details have been 
explained extensively in Ref. [ll|. Here, we just introduce 
some important propositions. 

In the E — H diagram, for each state {s^}, the energy E 
is represented by a straight line with slope — M = — Si 
since 



Proposition 1 For H < —h^i^^ (H > 
ground state is {si = —1} ({s^ = +1}). 



the 



with 



E{{s,},H)^Eo{{s^})-HM 



Eo{{si}) ^ - ^ J SiSj - ^ hi 

<i,j> i 



(4) 



(5) 



the energy axis intercept, i.e. the total energy of the 
configuration when H = 0. Consider a D-dimensional 
hypercubic lattice of size N — . Let /imax (^min) be 
the maximum (minimum) values of hi for a certain real- 
ization of the random fields. (Usually, random fields are 
chosen from a Gaussian distribution with mean and 
standard deviation R. R is often called the disorder pa- 
rameter.) Most of the following simple propositions have 
been proven in Ref. Tl. 



Proposition 2 Let the spin configuration Ci (C2) be the 
ground state for H — Hi (H — El'i). They have mag- 
netization Ml and M2, respectively. If Ci ^ C2 and 
H2 > Hi, then M2 > Mi. 

Thus, when sweeping the external field from H = —00 
to -ff = 00, the magnetization M will increase monotoni- 
cally. A corollary of this proposition is that: IntheE — H 
diagram if the slopes of the lines corresponding to the 
ground states Ci and C2 are different, i.e. Mi ^ M2, and 
without loss of generality we can assume Mi < M2, then 
the lines intersect at a field H* such that Hi < H* < H2 . 
This field H*{Ci, C2) is defined as the crossing field be- 
tween Ci and C2. According to the definition, one has 
Eq{Ci) - H*Mi = Eo{C2) - H*M2, so 



H*[Ci,C2) 



Eo{C2) — Eq{Ci 
M2 - Ml 



(6) 



For example, we can calculate the crossing field between 
the two simplest ground states: Ci = {si = —1} with 
Ml = -N and C2 = {s^_ = +1} with M2 = N. We have 
H* = —i/N'^^hi = —hi. For a zero-mean distribution 
of the random fields, we should have H* = 0. 

Proposition 3 Let the spin configuration Ci (C2) be the 
ground state for H = Hi {H = H2). Ci ^ C2, H2 > Hi 
and the crossing field between Ci and C2 is H* . If there is 
no configuration C such that E{C,H*) < E{Ci, H*) = 
E{C2, H*) then: (i) Ci is the ground state at least for 
the field range [Hi,H*) and (ii) C2 is the ground state 
at least for the field range {H* , H2]- 

This is the most important proposition. Its power 
comes from the fact that it can be applied iteratively. 
And finally we get the M{H) curve with all the ground 
states. See Fig. [1] 

Proposition 4 // the ground state is non-degenerate, 
then there can not he more than one avalanche connecting 
two nearest ground states in the E ~ H diagram. 

The proof of this interesting proposition is shown in 
the Appendix. This proposition says that if the GS is 
non-degenerate, when we increase the external field H 
adiabatically slowly, we can trigger just one avalanche at 
a time. 



III. PREPARATIONS OF THE PROOF 

In this section, we will work out the total energy change 
of the spin configuration due to multiple spin fiips and 
external field varying. The spin configuration is not nec- 
essarily the ground state. 

First, let's consider the simplest case of a single spin 
fiip. Suppose only one spin (s^) flips during the evolu- 
tion of configuration C at to configuration C at H' , 
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FIG. 1: An illustration of the algorithm to calculate the equilibrium M-H curve. Calculate the energies Ei and E2 of the two 
simplest states Ci = {si = —1} and C2 ~ {si = +1}, respectively, as a function of H. According to Proposition. [T] Ci (or C2) 
would be the ground state for H < — /imax (or H > — /imin). Calculate the crossing field H* {Ci,C2) where Ei = £2. Check 
whether there is a GS at H* which is different from Ci and C2. If no, the algorithm ends. If yes, denote the GS as C, calculate 
the crossing field H*{C, Ci) and H*{C, C2). At the new crossing fields, check whether there is a GS which is different from the 
two intersected states. The algorithm will not end until all the crossing fields have been checked. An example of the calculated 
equilibrium M-H curve is shown in Fig. (2] 
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H 

FIG. 2: The equilibrium M-H curve (the ground state evo- 
lution) for the Gaussian RFIM with _D = 3, L = 32 and 
R = 2.837. Here R is the standard deviation of the Gaus- 
sian random field distribution. The inset shows a detail of 
the M{H) curve near H = 0, where magnetization jumps are 
clearly seen. 



with AH ^ H' - H and AM = M' - M. Define Ui (or 
n^) to be the number of the i-th spin's nearest neighbors 
that point in the same direction as the spin in configu- 
ration C (or C"). We call these spins the same-direction 
nearest neighbors (SDNN) of the i-th spin. Note that 
Ui = 0, 1, 2, ...Z with Z = 2D the coordination number 
of the _D-dimensional hypercubic lattice. 



It is easy to get the bond energy change AJfjii — D). 
And the total energy change due to the single spin-flip 
and the varying external field is given by 

MH,AH) = MH)-AHM' (7) 

Here we have defined 

MH) = f,{H, 0) = 4J(n, - D) - [h, + H)As, (8) 

which is the energy change due to spin i flipping for the 
configuration C just at the field iJ, i.e. Asi = ±2 with 
AH = 0. It is easy to check that 

k±{H) = /.,± (0) ±2H = ±2hf (9) 

with '±' represents Si = ±1 and As; = =f2 accordingly. 

Second, we consider two spin flips. Suppose two differ- 
ent spins {si and Sj) flip during the evolution of configu- 
ration C at H to configuration C" at H' . There arc two 
subcases. 

(1) Si and Sj are not next to each other. The energy 
change is 

/,,,(H, AH) = MH) + fj{H) - AHM' (10) 

(2) Si and Sj are next to each other. The energy change 

is 

f(,^,){H,AH) = U{H)+f,{H)-U{s,-s,)-AHM' (11) 
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Note that the term —AJ{si ■ sj) is just due to the fact 
that the energy of the i — j bond will not change during 
the flip. 

Finally, let's consider the general case. See Fig. [3I^a). 
There are many spin flips during the evolution of config- 
uration C at field H to C at field H'. It is easy to check 
that the total energy change is given by 

AE{H,AH) = [f,{H) + fjiH) + ■ ■ ■] 
-4 J (sj • + • • • ) 
-AH{AI + As, + A.sj + • • • ) (f 2) 

On the RHS, the first term includes all the flipping spins. 
The second term includes all the nearest-neighbor inter- 
actions among those flipping spins. The last term is due 
to the varying external field. In particular, if all the flip- 
ping spins flip at the same H and they are connected 
to each other and have the same spin value —1 (or +1) 
before the flip, then this collective spin flip is called an 
avalanche (or a reverse avalanche). 

Denote the energy change due to an avalanche as 
fA„{H,AH), we have 

fA^ {H, AH) = [U{H) + f,{H) + ■■■]- 4J7Vb(A„) 
-AiJ(M + 2Sa) 
= fAAH)-AH{M + 2So,) (13) 

with iVb(Ac() defined as the number of interacting bonds 
in Aa, Sa the size of the avalanche and fA^{H) the en- 
ergy change due to the avalanche when AH = 0. Simi- 
larly, for the reverse avalanche, we have 

(i/, AH) = [MH) + MH) + ■■■]- 4J7Vb(A^) 
-AH{M ~ 2S'p) 
^ fA^^^{H)-AH{M-2Sl) (14) 

Due to Eq. [9l we have 

fAjH) = fAA0)-2S,,H (15) 
/^.(i/) - fA'^iO) + 2S'pH (16) 

Now we can rewrite Eq. [1^ in terms of /a„ and Ja'^ ■ 
The total energy change due to avalanches and reverse 
avalanches is given by 



IV. PROOF OF THE NO-PASSING RULE 

Now we are ready for the proof of the no-passing rule. 
Let the spin configuration Ci (C2) be the ground state for 
H = Hi {H = H2). H2> Hi. Suppose Ci and C2 are 
connected with multiple avalanches: A\,A2- ■ - An with 
sizes 5*1, S'2 • • • Sn and reverse avalanches • • • AJ!^, 

with sizes S\.,S\^--- respectively. To compensate 
these reverse avalanches (so as to make sure M is mono- 
tonically increasing, see Proposition. [2]), we must have 
Sa = ELi So. > Sa^ = See Fig.^a). 

The idea is that if C2 is the GS at field H2, then it 
should have lower energy than any other spin configura- 
tion at H2. But we will prove this is NOT true. Just 
consider another spin configuration C. The only differ- 
ence between C2 and C is that C is evolved from Ci 
without any reverse avalanches. See Fig. [3]Jb). We now 
try to prove that E{C, H2) < E{C2, H2), so C2 can NOT 
be the GS at H2. But this is equivalent to proving that 
AE < AE. Here, 

AE = E{C2,H2)-E{Ci,Hi) 

= FAr{Hi) + Fa{Hi) + 4 J iVb(^, A') 

- AH{M -2SAr + 2Sa). (18) 

On the other hand, 

A^ = E{C,H2)-E{Ci,Hi) 

= Fa{Hi)-AH{M + 2Sa) (19) 

Therefore, 

AE-AE = Fa'{Hi)+UNi,{A,A')+2Sa'AH > (20) 

Here we have used the fact that Ci is the ground state 
for H = Hi such that any kinds of spin flip will increase 
the energy: fA'^{Hi) > =^ Fa^{Hi) > 0. Also, for the 
ferromagnetic RFIM, J > 0. Since each term is positive, 
so the sum is positive, i.e. AE > AE or E > E. Actu- 
ally, for any state C2 which evolved from Ci with reverse 
avalanches, we can find a corresponding state C which 
evolved from Ci without any reverse avalanches that has 
lower energy than C2 at field H2. So reverse spin flips 
are impossible for ground state evolution when increas- 
ing external fleld. Generally, flipped spins can never flip 
back when we sweep the external field monotonically. 



AE{H,AH) = FA{H)+FAr{H)+ AJNi,{A,A') 

- AH{M + 2Sa - 2Sa') (17) 

with notations Fa{H) = Y^c^JaJH), Fa^{H) = 
EpfA'.iH), Sa ^ EcSc, Sa^- ^ Y.pSl. Here 
A'^b(^, ^'^) denotes the number of interacting bonds be- 
tween avalanches and reverse avalanches. For example, 
in Fig. Ha), iVb(A^O = 3. 



V. APPLICATION 

The straightforward application of the no-passing rule 
is very useful to accelerate the calculation of the ground 
states when varying the external field. Suppose the GS 
Ci at field Hi has already been obtained, and we want 
to calculate the GS C2 at field H2 with H2 > Hi. Ac- 
cording to the no-passing rule, the UP spins in Ci will 
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FIG. 3: States evolved from Ci, with only the change of the spin configuration, i.e. avalanches and reverse avalanches are 
explicitly shown, (a) State C2: evolved from state Ci with both avalanches and reverse avalanches. (Black dot) spins flip UP, 
forming avalanches {Ai, A2, A3 and A4). (White dot) spins flip DOWN (reverse flip), forming reverse avalanches {A{, A2). 
Note that there are three interacting bonds between avalanche A2 and reverse avalanche A2- (b) State C evolved from state 
Ci without reverse avalanches. 



stay UP in C2, i.e. those spins are frozen, so we needn't 
consider them in the ground state analysis. We just need 
to consider the DOWN spins in Ci. The only cost is 
that we have to deal with the frozen UP spins as com- 
plicated fixed boundary conditions for the DOWN spins. 
At first sight, one may think that only when the den- 
sity of the frozen spins is big enough can we make the 
GS calculation faster. But how big is enough? To opti- 
mize our calculation, we consider the running time dif- 
ference (At) between the two methods: (I) without using 
the earlier solution Ci ; (II) using the earlier solution Ci . 
For both methods, ground states are found by using the 
push-relabel algorithm. The numerical experiments are 
conducted on a desktop with 2.80 GHz CPU and 2GB 
Memory. And we tune the UP-spin density riup (Down- 
spin density ridown) by varying Hi . The result is shown in 
Fig.m It is found that for H2 > Hi, as long as n^p > 0.07 
in GS Ci, method II will be faster than I. Symmetrically, 
for H2 < Hi, as long as ndown ^ 0.07 in GS Ci, method 
II will be faster than I. This suggests it is not necessary to 
have an extremely large portion of frozen spins to use the 
earlier solution. Freezing a tiny part of spins will accel- 
erate the GS calculation already. Furthermore, for larger 
and larger density of the frozen spins, using the earlier 
solution will save more and more running time. (Keep 
in mind that for RFIM, the running time of the push- 
relabel algorithm scales as 0{N'^^'^).) Consequently, the 
calculation of the whole M{H) curve will be sped up 
dramatically. 



VI. DISCUSSIONS 

Throughout our proof of the no-passing rule, we don't 
assume that the ground state is unique. In other words, 
the no-passing rule is correct even when the ground state 
is degenerate. For example, this happens for the RFIM 
when the random fields are chosen from a bimodal dis- 
tribution [1^ . 

In the proof we explicitly use the fact that the nearest 
neighbor interaction should be ferromagnetic (J > 0). 
This means any antiferromagnetic interactions will de- 
stroy the no-passing rule. Thus, for other random mag- 
net models, if Jij could be negative, such as the random- 
bond Ising model (RBIM) with negative Jij or the spin 
glasses, the rule will be violated. 

Finally, we conjecture that for elastic manifolds in 
random media, there could be a similar equilibrium no- 
passing rule at zero temperature, provided that the elas- 
tic potential is convex and partial ordering of the config- 
urations can be clearly defined. 
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FIG. 4: Running time difference (At) between metliods witfi 
and witliout using tlie earlier solution, i.e the GS Ci at field 
Hi, to calculate the GS C2 at field H2- The time difference At 
(given in seconds) is plotted against the UP-spin (or DOWN- 
spin) density of the GS Ci. At < means using the earlier 
solution will save the running time. Calculations are done 
for 3D Gaussian RFIM (with disorder parameter R = 2.27) 
for different system sizes. (Top) H2 > Hi. UP spins in Ci at 
field Hi will stay UP at field H2. (Bottom) H2 < Hi. DOWN 
spins in Ci at field Hi will stay DOWN at field H2. 



APPENDIX A: A SINGLE AVALANCHE 
CONNECTS TWO NEAREST GROUND STATES 

Here, we show the proof of Proposition. 0) 



Proof: Suppose when the field is increased from Hi 
to H2, the GS Ci evolves to the nearest GS C2 with two 
avalanches {Ai and A2 with size 6*1 > 1 and 6*2 > 1, 
respectively) . 

The crossing field is given by 



Eq{C2) — Eq{Ci 



M2 



2iSi^S2) 

The last line is due to Eq[T7l We can choose a trial 
state C which is evolved from Ci with only avalanche Ai 
occurring. We want the following relation to hold 



E{C,H*)<E{CuH*) = E{C2.H*) (A2) 
which is equivalent to 

E{C,H*)-E{CuH*) = fAdH*) 

< (A3) 

Plugging Eg. lAll in it, we just need to prove 



(A4) 



If we do have 5*2/711 (0) < SiJa^ (0), then we choose the 
trial state C, which has lower energy than Ci and C2 at 
the field H* ; if 6*2 /ai (0) > SifA2 (0), then we can choose 
another trial state C", which is evolved from Ci through 
only avalanche A2 and has lower energy than both Ci and 
C2 at H* . In both cases, we have shown that C2 can not 
be the nearest GS at H2 for the GS Ci at Hi, if Ci evolves 
to C2 with two avalanches. If S2fAA^) = "51/^2(0), it is 
easy to show that £;(Ci, i?*) = E{C,H*) = E{C',H*) = 
E{C2, H*), then there will be degenerate ground states at 
H* , which is in contradiction to the hypothesis of Propo- 
sition, m 
Q.E.D 
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